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ABSTRACT 

Submillimeter Array observations of Orion-KL at ~ 1" resolution in the 440 /im/690 GHz band reveal 
new insights about the continuum and line emission of the region. The 440 /Ltm continuum flux density 
measurement from source I allows us to differentiate among the various proposed physical models: Source 
/ can be well modeled by a "normal" protostellar SED consisting of a proton-electron free-free emission 
component at low frequencies and a strong dust component in the submillimeter bands. Furthermore, 
we find that the protostellar object SMA1 is clearly distinct from the hot core. The non-detection of 
SMA1 at cm and infrared wavelengths suggests that it may be one of the youngest sources in the entire 
Orion-KL region. The molecular line maps show emission mainly from the sources I, SMA1 and the hot 
core peak position. An analysis of the CH 3 CN(37_k: — 36^) iGladder (K — 0...3) indicates a warm gas 
component of the order 600 ± 200 K. In addition, we detect a large fraction (~ 58%) of unidentified lines 
and discuss the difficulties of line identifications at these frequencies. 

Subject headings: techniques: interfcrometric — stars: formation — ISM: individual (Orion-KL) — 
ISM: molecules — ISM: lines and bands — submillimeter 



1. INTRODUCTION 

The 440^m/690GHz band is almost entirely unexplored 
at high spatial resolution. In February 2005, the Submil- 
limeter Array (SMA 3 ) performed the first imaging cam- 
paign in that band with 6 antennas in operation, thus line 
and continuum imaging at 1" resolution was possible at 
frequencies around 690 GHz. One of the early targets of 
this campaign was the Orion-KL region. 

At a distance of ~ 450 pc, Orion-KL is the closest 
and most studied region of massive star formation in our 
Galaxy. Furthermore, the Orion-KL hot core is known 
to be particularly rich in molecular line emission (e.g., 
Wright et al. 1996; Schilke et al. 2001; Beuther et al. 2005). 
The hot core as traced by (sub)millimeter dust continuum 
emission is a chain of clumps offset by ~ 1" from the radio 
source /. The hot core is rich in nitrogen-bearing molecules 
like CH3CN and NH3 and has temperatures of the order 
of a few hundred Kelvin (e.g., Wilner et al. 1994; Wilson 
et al. 2000; Beuther et al. 2005). In contrast, the so-called 
compact ridge approximately 5" to the south is partic- 
ularly rich in CH3OH emission and has lower tempera- 
tures of the order 100 K (e.g., Wright et al. 1996; Liu et al. 
2002; Beuther et al. 2005). In addition, the region harbors 
a complex cluster of infrared sources studied from near- 
to mid-infrared wavelengths (e.g., Dougados et al. 1993; 
Gezari et al. 1998; Greenhill et al. 2004; Shuping et al. 
2004). At least two outflows are driven from the region, 
one high-velocity outflow in the south-east north-west di- 
rection (e.g., Allen & Burton 1993; Chernin & Wright 
1996; Schultz et al. 1999), and one lower- velocity outflow in 
the north-east south-west direction (e.g., Genzel & Stutzki 



1989; Blake et al. 1996; Chrysostomou et al. 1997; Stolovy 
et al. 1998; Beuther et al. 2005). The driving source(s) of 
the outflows are uncertain. Initial claims that it might be 
IRc2 are outdated now, and possible culprits are the radio 
source / and/or the infrared source n, also known as radio 
source L (Menten & Reid 1995). 

In 2004, we observed Orion-KL at 1" resolution with 
the SMA in the 865/mi/348GHz band (Beuther et al. 

2004, 2005). The molecular line data show a rich line 
forest resolving different spatial structures for the vari- 
ous species. The SiO(8-7) observations convincingly show 
that the north-east south-west outflow emission originates 
at source /. Temperature estimates based on a CH3OH 
multi-line fit result in values as high as 350 K for the hot 
core. In the 865 fj,m submillimeter continuum data, source 
/ was resolved from the hot core and source n was detected 
as well. The spectral energy distribution of source / from 
8 to 348 GHz still allowed ambiguous interpretations of 
its physical nature (sec also §3.1.3). Furthermore, a new 
continuum peak, SMAl, between sources / and n was de- 
tected; however, the nature of SMAl was uncertain, as it 
could have been either an independent protostellar core or 
part of the extended hot core. In this paper, our new ob- 
servations in the 440 /im band clarify the nature of SMAl 
and source /. 

2. OBSERVATIONS 

We observed Orion-KL with the SMA on February 19th, 

2005, in the 440/im/690GHz band. Six antennas were 
equipped with 690 GHz receivers, and the covered base- 
lines ranged between 16 and 68 m. The weather condi- 
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tions were excellent with a zenith opacity, - measured 
with the NRAO tipping radiometer located at the Caltech 
Submillimeter Observatory - r(230GHz) between 0.03 and 
0.04 throughout the night. This corresponds to zenith 
opacities at 690 GHz between 0.6 and 0.8 (t(690GHz) ~ 
20 x (t(230GHz) - 0.01), Masson 1994). The phase center 
was the nominal position of source / as given by Plambeck 
et al. (1995): R.A. [J2000] 5 h 35 m U.50 s and Dec. [J2000] 
-5°22'30".45. The FWHM of the primary beam of the 
SMA at these frequencies is 18", the aperture efficiency 
is ~ 0.4. For more details on the array and its capabili- 
ties see Ho et al. (2004), and recent updates are given at 
http://smal.sma.hawaii.edu/specs.html. 

The receivers operated in a double-sideband mode with 
an IF band of 4-6 GHz so that the upper and lower side- 
bands were separated by 10 GHz. The correlator had 
a bandwidth of 1.968 GHz and the channel spacing was 
0.8125 MHz. The covered frequency ranges were 679.78 to 
681.75GHz and 689.78 to 691.75GHz (Fig. 1). Bandpass 
calibration was done with observations of Callisto (angu- 
lar diameter ~ 1.4"). Since no quasar is strong enough 
in that band for calibration purposes, we calibrated phase 
and amplitude via frequent observations of Titan (angu- 
lar diameter <~ 0.9"). The distances of Callisto and Ti- 
tan from Jupiter and Saturn were > 230" and ~ 190", 
respectively. The contamination of the Callisto and Ti- 
tan fluxes by their primaries is negligible due to mainly 
three reasons: (a) the primary beam pattern severely at- 
tenuates these planets since they are much more than 10 
primary beams away, (b) the large size of the planets (40" 
for Jupiter, 20" for Saturn) implies that a large fraction of 
their fluxes is resolved out, and (c) Due to the motion of 
the moons over the period of observations any small con- 
tamination will not be added coherently further reducing 
the contamination issue. 

The large angle between Titan and Orion-KL (~ 40°) 
introduced larger than normal uncertainties in the phase 
and amplitude transfer. The flux calibration was per- 
formed with the Titan observations as well, for which a 
well-calibrated model has been created and used exten- 
sively at mm wavelengths (Gurwell & Muhleman 1995, 
2000; Gurwell 2004). In the current observations, the up- 
per sideband (USB) included the CO (6-5) line, which is 
quite broad as determined by using both observations of 
lower CO transitions and modeling of the expected line 
shape, assuming that CO has a uniform 50 ppm abun- 
dance in the Titan atmosphere (Gurwell 2004). Modeling 
the flux of the continuum LSB and continuum+CO(6-5) 
USB emission, the expected average Titan fluxes at the 
given dates were 12.0 and 22.3 Jy for the LSB and USB, 
respectively. Based on flux measurements of the two main 
continuum sources in Orion-KL using only the part of the 
spectra with no detectable line emission (§3.1), we esti- 
mate the relative flux accuracy between the upper and 
lower sideband to be better than 10%. The absolute flux 
calibration was estimated to be accurate within ~ 25%. 

Since Orion-KL is close to the celestial equator, and we 
have for most baselines only source data observed in an 
hour angle range from —2.6 to +1.2 hours (the data of 
only 3 baselines were usable over a longer range from —2.6 
to +4.5 hours), imaging of the data proved to be difficult. 
Figure 2 shows the given uv-coverage and the resulting 



dirty beam. In spite of the large-scale emission present 
in Orion-KL, we had to employ clean-boxes covering the 
central region around sources I, SMA1 and the main hot 
core peak to derive reasonable images because of the high 
side-lobe levels. Measured double-sideband system tem- 
peratures corrected to the top of the atmosphere were be- 
tween 900 and 3500 K, depending on the elevation of the 
source. Our sensitivity was limited by the strong side-lobes 
of the strongest emission peaks, and thus varied between 
the continuum and the line maps of different molecules 
and molecular transitions. This limitation was mainly due 
to the incomplete sampling of short uv-spacings and the 
presence of extended structures. The achieved la rms 
of the 440 /im continuum image, produced by averaging 
the apparently line-free part of the USB spectrum (see 
§3.1), was ~ 700mJy/beam - well above the expected 
lcr of ~ 41mJy/beam for the given on-source time, the 
used bandpass and an average r(230GHz) = 0.037 (cor- 
responding to 0.61mm precipitable water vapor). This 
larger rms was mainly due to the strong side-lobes and 
inadequate cleaning. The theoretical la rms per lkms^ 1 
channel was ~ 775mJy/beam, whereas the measured la 
rms in lkms -1 channel images was ~ 2.6 Jy/beam, again, 
because of the strong side-lobes and inadequate cleaning. 
The la rms for the velocity-integrated molecular line maps 
(the velocity ranges for the integrations were chosen for 
each line separately depending on the line-widths and in- 
tensities) ranged between 1.1 and 1.6 Jy/beam. The syn- 
thesized beams were 1.4" x 0.9" using an intermediate 
weighting between natural and uniform ("robust" value 
in MIRIAD). We calibrated the data within the IDL 
superset MIR developed for the Owens Valley Radio Ob- 
servatory and adapted for the SMA (Scoville et al. 1993, 
see also the MIR cookbook by Charlie Qi at http://cfa- 
www.harvard.edu/~cqi/micook.html). The imaging was 
performed in MIRIAD (Sault et al. 1995). 

3. RESULTS AND ANALYSIS 

3.1. 440iim continuum emission 

The spectra presented in Figure 1 exhibit many molec- 
ular lines, therefore constructing a continuum image had 
to be done with great care. The continuum image pre- 
sented in Fig. 3 (right panel) was constructed by averag- 
ing the part of the USB spectrum with no detectable line 
emission as shown in Fig. 1. The corresponding appar- 
ently line-free part of the LSB is considerably smaller and 
thus the comparable pseudo-continuum image noisier. Al- 
though we therefore omit the LSB data for the scientific 
interpretation, we point out that the peak fluxes in the 
LSB continuum image vary by less than 10% compared to 
the USB image. This shows the good relative flux calibra- 
tion between both sidebands. 

Weak lines below the detection limit may contribute to 
the continuum emission. To estimate the contamination 
from such weak lines, we produced an additional pseudo- 
continuum dataset of the lower sideband data including 
all spectral lines. We used the LSB data for this purpose 
because its line-contamination is more severe. Fitting two 
point-sources in the uv-domain (a) to the whole LSB con- 
tinuum data and (b) to the apparently line-free LSB data 
(for details see §3.1.1), the measured fluxes from (a) were 
only ~ 15% higher for source / and SMA1 than those from 
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(b). In our data, possible weak lines unidentified due to 
noise are on average less than 10% of the average peak-flux 
of the strong lines (Fig. 1). Since the strong lines already 
cover about half of the given bandpass (Fig. 1), the contri- 
bution from possible weak lines below the detection limit 
amounts to only about 1.5%. This is negligible compared 
to the absolute flux uncertainties. 

The derived 440 /im continuum emission map is shown 
in the right panel of Figure 3. To make a proper com- 
parison with the 865 /im emission, we also show the orig- 
inal image at that wavelength (Beuther et al. 2004), as 
well as an image of the 865 observations using only 
the 30 to 160 kA data range, which is the range of base- 
lines sampled by our 440 /im data. In spite of strong 
side-lobes and the imaging problems described in §2, the 
440 /im map shows two distinct point-like sources associ- 
ated with source / and SMA1. Comparing this image with 
the 30 to 160 kA image at 865 /im, we find at the longer 
wavelength similar structures for source / and SMAl plus 
emission from the hot core. Additional features like source 
n are not present in the 440 /im continuum image, likely 
because of the poor imaging capabilities and thus strong 
side- lobes of the observations (sec §2). The main difference 
between the 440 /im and the 865 /im continuum images is 
the strong presence of the hot core in the longer wave- 
length image, whereas we find only a small extension to 
the east at 440 /im. Imaging of the spectral line data in the 
440 /im band (§3.2) shows that, in spite of the mentioned 
imaging limitations (§2), we can spatially image and sep- 
arate line emission from source / and the hot core. Hence, 
while it remains puzzling that the typically strong hot core 
is only barely detectable in our 440 /im image, it might be 
a real feature of the data indicating a peculiar low spectral 
index a of the hot core emission (for a detailed discussion 
see §3.1.5). 

3.1.1. Flux measurements 

The 440 /im continuum emission in Fig. 3 gives the im- 
pression that the point sources / and SMAl are embedded 
in a larger-scale ridge emission. While Orion-KL is known 
to be part of the large-scale Orion molecular ridge (e.g., 
Lis et al. 1998), the larger-scale emission in our image is 
to some degree also an artifact of the data reduction, since 
we had to use clean boxes and thus shifted emission from 
the entire primary beam to the central region. Hence, flux 
measurements toward source / and SMAl from the images 
will give values that are too high. To avoid these problems, 
we fitted the data in the uv-domain assuming: 

(a) two point sources at the positions of source / and 
SMAl; 

(b) two point sources at the positions of source /and SMAl 
+ an elliptical Gaussian; 

(c) three point sources at the positions of source /, SMAl, 
and the hot core peak position; 

(d) three point sources at the positions of source /, SMAl, 
and the hot core peak position + an elliptical Gaussian. 
Afterwards, we imaged the models using exactly the same 
uv-coverage and imaging parameters employed previously 
for the real data. Figure 4 shows a direct comparison be- 
tween the image of the original data and three model fits. 
We find that the models with additional larger-scale ellip- 
tical Gaussian component reproduce the data rather well. 



Although the fitted fluxes for source / and SMAl vary by 
less than 15% (within the calibration uncertainty) with or 
without the third point source at the position of the hot 
core, the three-point-sources+Gaussian model image most 
resembles the original data image. Therefore, we use this 
model for the following analysis. 

These simplified models allow us to bracket the fluxes 
from sources /, SMAl and the hot core. The three-point- 
source-only model likely overestimates the intrinsic fluxes 
because it adds some of the underlying large-scale emission 
to the point sources. Contrary to this, the three-point- 
source+Gaussian model underestimates the point-source 
contributions and thus the point-source fluxes. The mea- 
sured values for both uv-fitting approaches are: 



S(I) 

S*(SMA1) 

S*(HotCore(HC)) 

Gaussian 



3-point 3-point+Gaussian 
7.9 Jy 4.6 Jy 

7.1 Jy 3.5 Jy 

3.1 Jy 2.4 Jy 

178 Jy (size 6.8" x 2. 



Including the 25% calibration uncertainty, the 440 /<m 
fluxes measurements with conservative estimates of the er- 
ror budget are: 

3.5 < S(I) < 9.9 => S(I) ~ 6.7 ± 3.2 Jy 

2.6 < S(SMA1) < 8.9 S(SMA1) - 5.75 ± 3.15 Jy 

1.8 < S(HC) < 3.9 S(HC) ~ 2.85 ± 1.05 Jy 

The errors are not to be understood as la or 3a values, 
but they rather give the extreme values of what the fluxes 
of the three sources could be. As outlined above, the line 
contamination of our continuum dataset is very low. 

3.1.2. The large-scale structure 

Obviously, from our data we cannot constrain well the 
large-scale flux distribution. However, the fitted 178 Jy 
are a lower limit to the large-scale flux since we have no 
data below 28 kA, and Johnstone & Bally (1999) measured 
a 450 /urn single-dish peak flux of 490 Jy with SCUBA at 
7.5". Furthermore, we investigated the influence of a large 
source with a flux of a few hundred Jy on our observa- 
tions. Modeling different source sizes and observing them 
with the given uv-coverage, the size-scale of the large-scale 
emission has to exceed 11". Larger sources are filtered out 
by the observations and do not distort the observed con- 
tinuum image significantly, whereas smaller sources with 
large flux would completely dominate the observations. 
The scale of 11" corresponds approximately to the size- 
scales theoretically traceable by the shortest baselines of 
28 kA (corresponding to approximately 9" scales). 

3.1.3. Source I 

Prior to these observations, the SED of source / could 
be fitted well by two very different models, either as the 
result of optically thick proton-electron free-free emission 
up to 100 GHz plus dust emission that accounts for most 
of the submm flux, or H~ free-free emission that gives rise 
to a power-law spectrum with a power-law index of ~ 1.6 
(Beuther et al. 2004). Since the two models vary strongly 
at higher frequencies, an accurate 440 /im continuum flux 
measurement was of utmost importance in discriminating 
between these models. 
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Figure 5 shows the complete SED of source / from 8 
to 690 GHz. The measured 690 GHz flux fits the model 
of a proton-electron free-free plus dust emission spectrum 
well, but not the model of the power-law spectrum ex- 
pected for H~ free-free emission. Therefore, the 440 fim 
continuum data can discriminate between the two models, 
and the SED of source / turns out to be a rather typical 
SED for deeply embedded protostellar sources. The im- 
portant result here is that the turn-over frequency (i.e., 
the transition from optically thick to optically thin emis- 
sion) is much higher in source / than in other protostars 
where it is typically in the range of 10 to 40 GHz. 

3.1.4. SMAl 

The observations show that SMA1 is not just an exten- 
sion of the hot core, but is an independent protostellar 
source. Based on PdBI HC3N(10-9) observations of vi- 
brational excited lines, de Vicente et al. (2002) predicted 
an additional massive protostellar source in the region ap- 
proximately 1.8" south of source /, which is the location 
of SMAl. Only two flux measurements exist for SMAl: 
at 865 /xm (S ~ 0.36 ± 0.07 Jy, Beuther et al. 2004) and at 
440 A*m (S — 5.75 ± 3.15 Jy). The spectral index between 
these two bands is S oc v a with a ~ 4.0^°'4- In spite of the 
broad range of potential a, the data are consistent with a 
typical dust opacity index (3 — a — 2 in star-forming re- 
gions of 2. The non-detection of SMAl at cm and infrared 
wavelengths suggests that it may be one of the youngest 
sources of the evolving cluster. 

3.1.5. The Hot Core 

The hot core flux values we derived above are probably 
more uncertain than those from source / and SMAl be- 
cause the hot core is only weakly discernible from the rest 
of the image in the original data as well as the model. Fur- 
thermore, the response of the continuum and line emission 
to the poor uv-sampling is not the same in interferometric 
studies, and we cannot entirely exclude that the weakness 
of the hot core in the 440 /im data is caused by this poor 
uv-sampling and the resulting imaging problems. Keep- 
ing in mind these uncertainties, we nevertheless can use 
the existing data to estimate a spectral index between the 
previous 865 /xm measurement and this new data. The flux 
measurements in the two bands toward the hot core peak 
position (S"(865/zm) ~ 0.54±0.11Jy, Beuther et al. 2004, 
and 5(440/im) ~ 2.85 ± 1.05 Jy) result in a range of po- 
tential spectral indices S oc v a with a ~ 2.4+g'g. These 
values are considerably lower than for source / and SMAl. 

Single-dish spectral index studies toward the Orion A 
molecular cloud qualitatively find a somewhat lower spec- 
tral index toward Orion-KL than toward most of the rest 
of the Orion A molecular cloud (Lis et al. 1998; Johnstone 
& Bally 1999), but the two studies do not agree quantita- 
tively. Lis et al. (1998) find a dust opacity index [5 ~ 1.8 
corresponding in the Raylcigh- Jeans limit to a spectral 
index a ~ 3.8, whereas Johnstone & Bally (1999) find a 
spectral index a ~ 2.2. Nevertheless, considering the large 
number of potential uncertainties (e.g., calibration, differ- 
ent beam sizes, line contamination, etc.), the qualitatively 
lower spectral index a toward Orion-KL compared with its 
surroundings appears a reliable result from the single-dish 
studies. 



In contrast to source / and SMAl, where the spectral 
indices are both consistent with typical values of a around 
4 in star-forming regions, the lower spectral index toward 
the hot core is intriguing. Furthermore, lower-resolution 
observations show that on larger spatial scales the data are 
dominated by the apparently more extended hot core emis- 
sion (e.g., Plambeck et al. 1995; Blake et al. 1996; Beuther 
et al. 2004). Therefore, the different spectral indices indi- 
cate that the lower single-dish spectral index may be due 
to hot core emission distributed over larger spacial scales, 
whereas the protostellar sources / and SMAl exhibit spec- 
tral indices more typical of those from other star-forming 
regions. 

3.2. Spectral line emission 

We detected 24 spectral lines over the entire bandpass 
of 4 GHz (Fig. 1 & Table 1). This is a significantly lower 
line number than the ^150 lines over the same spectral 
width at 348 GHz. The lower sensitivity of the 690 GHz 
data compared with the 345 GHz observations accounts for 
this difference to some degree. Furthermore, at 690 GHz 
each line of a given Aw covers twice the frequency range 
than at 345 GHz, and thus strong lines dominate a broader 
part of the spectrum than at lower frequencies. Further- 
more, to first order the observing frequency is proportional 
to the average energy levels of observable lines in that 
frequency range, and the lower number of line detections 
at 690 GHz indicates that the bulk of the gas is at lower 
temperatures (see temperature discussion in §3.3). Nev- 
ertheless, the energy levels of the 690 GHz lines (Table 1) 
show that a very warm gas component has to be present 
in the region as well. Temperature estimates based on the 
CH 3 CN(37k — 36^) are derived below (§3.3). Compar- 
ing these spectra with the single-dish spectra toward the 
same region (Schilke et al. 2001; Harris et al. 1995), we 
find a large number of lines previously not detected in this 
region. Line identifications were done first using the line 
survey by Schilke et al. (2001), and then cross-checking 
the identifications with the molecular databases at JPL 
and CDMS (Poynter & Pickett 1985; Miiller et al. 2001). 
It has to be noted that the unambiguous identification of 
spectral lines becomes more difficult at higher frequencies, 
which is discussed in §3.4. 

Table 1 also shows the line-widths derived from Gaus- 
sian fits to the vector-averaged spectra on a short base- 
line of 25 m (Fig. 1). Since the measured line- widths are 
associated with gas from the whole emission in the pri- 
mary beam, they do not represent the line-width toward 
a selected position but rather an average over the whole 
region. The derived values for Av range between 4.8 and 
16.4kms _1 . We could not fit the CO(6-5) line due to 
the missing flux. While some of the broader line-widths 
may be due to the outflows in the region, this is prob- 
ably not the main reason for most of the lines since the 
spectra often resemble a Gaussian profile. It seems more 
likely that different velocity components as well as large in- 
trinsic line-widths toward each position due to turbulence 
and/or internal motion cause a reasonable fraction of the 
line-width. For a discussion of the CH3CN line-widths, see 
§3.3. 

Figure 6 presents the images of the velocity-integrated 
emission of a few of the identified molecular lines in 



SMA 440^m/690GHz observations of Orion-KL 



5 



the bandpass. As discussed in §2 and §3.1, imaging of 
the emission turns out to be difficult , as evidenced by 
the strong negative side-lobes in Fig. 6. Since our pri- 
mary beam is rather small (18") and the imaging gets 
worse towards the edge of the primary beam, we consider 
only emission at the centers of the maps to be reliable. 
All line images in Fig. 6 show considerable emission to- 
ward source 7, SMA1 and the hot core peak position. 
C 33 S(14-13) peaks toward SMAl with lower level emis- 
sion from the hot core and source 7. In contrast to this, 
the line images from CH 3 CN(37 2 -36 2 ), CH 3 OH(22 1 -21 2 ) 
and S0 2 (353,33 — 34 2 ,3 2 ) all show emission peaks from 
source 7, SMAl and the hot core peak position. The 
HNCO(31 2 , 2 g — 30 2j2 8) image exhibits emission peaks as- 
sociated with source 7 and SMAl, and only lower- level 
emission from the hot core position. A striking difference 
between this 440 /im line images and our previous study in 
the 865 /jm band (Bcuthcr et al. 2005) is that the higher 
excited lines in the 440 /im band show line emission from 
source I. In the 865 /im band, we found mainly SiO emis- 
sion from source 7, other species appear to avoid the source 
in the lower frequency band. This indicates the presence 
of high temperatures in the close environment of source I. 

We could not produce reliable emission maps for the 
CO (6-5) and H 13 CN(8-7) lines. While this is less of a sur- 
prise for CO, which is known to be extended and thus dif- 
ficult to image interferometrically, H 13 CN was expected to 
be far more compact, comparable to many of the other im- 
aged lines. However, comparing all lines in our frequency 
setup (excluding only the CO(6-5) line), the H 13 CN(8-7) 
line has one of the broadest lincwidths (Av ~ 16.2 km s -1 ), 
and thus may be more strongly affected by the molecu- 
lar outflows present than initially expected. This could 
explain our inability to properly image the H 13 CN(8-7) 
emission. 

3.3. CH 3 CN(37k — 36^)' temperature estimates and 
line- widths 

Since we observed the K = 0,1,2,3 lines of the 
CH 3 CN(37k - 36 K ) lines (the K = 0,1 lines are blended), 
we can use this TT-series for temperature estimates of the 
warm gas component in the regions. Figures 6 & 7 show 
the integrated CH 3 CN(37 2 — 36 2 ) emission map and spec- 
tra of all TT-components toward the CH3CN peak positions 
associated with source 7, SMAl, and the hot core. 

Using the XCLASS superset to the CLASS software 
(Peter Schilkc, priv. comm.), we produced model spectra 
in local thermodynamic equilibrium (Fig. 7). This soft- 
ware package uses the line catalogs from JPL and CDMS 
(Poynter & Pickett 1985; Miiller et al. 2001). Because 
the K = 3 component has double the statistical weight 
than the K = 0, 1, 2 lines, one expects in the optical thin 
case larger line intensities for the K = 3 transition. The 
optically thin model spectra reproduce the relative line 
intensities well toward SMAl and the hot core. Toward 
source 7, the modeled K = 3 intensity is a bit larger 
than observed indicating increasing optical depth there. 
The upper energy levels E u of the various K-ladder lines 
range only between 621 and 685 K (Table 1), and the data 
are not very sensitive to the temperature. We find rea- 
sonable model spectra in the 600 K regime with a large 
potential range of temperatures of ±200 K. The CH3CN 



model column densities we used are of the order a few 
times 10 15 cm~ 2 , and the line- width adopted for the mod- 
cling was 6.5 km s -1 . In spite of the large errorbars, the 
temperatures traced by these high-energy lines are, as ex- 
pected, higher than the highest temperatures of 350 K pre- 
viously derived from our CH3OH multi-line analysis in the 
865 fim band (Beuther et al. 2005). We attribute this tem- 
perature difference mainly to the various gas components 
at different temperatures and densities, and to the differ- 
ent excitation energies of the spectral lines in the 690 and 
337 GHz bands. 

Furthermore, we can compare the modeled 
CH 3 CN(37k - 36^) line-widths of the order e^kms" 1 
with the previous SMA observations of the CH3CN(19 8 — 
18 8 ) line at 348 GHz (Beuther et al. 2005). The measured 
CH 3 CN(19 8 -18 8 ) from SMAl and the hot core are 4.6 and 
4.1kms _1 , respectively. Beuther et al. (2005) found only 
very weak CH3CN(19 8 — 18 8 ) toward source 7. Although 
the excitation temperatures E u of the CH3CN(37 2 — 36 2 ) 
and the CH 3 CN(19 8 — 18 8 ) lines are approximately the 
same (649 and 624 K, respectively), it is interesting that 
the line-width toward the CH 3 CN(37 2 — 36 2 ) is larger. 
Since the critical densities between the CH3CN(37k~ 36k) 
and CH 3 CN(19k — 18k) series vary by about one order 
of magnitude (~ 10 8 and ~ 10 7 cm" 3 , respectively), this 
indicates that the higher density gas CH 3 CN(37k — 36k) 
lines arc likely subject to significantly more turbulent mo- 
tions. As already pointed out by Sutton et al. (1986), 
radiative excitation is unlikely to significantly affect the 
populations of the vibrational ground state transitions, 
even for high J-levels. 

3.4. Unidentified lines (UL) 

Fourteen out of the 24 detected spectral lines remain 
unidentified or only tentatively identified. Such a high 
percentage of unidentified lines (~ 58%) has rarely been 
reported, to our knowledge only recently toward SgrB2 
(Friedel et al. 2004). The previous single-dish line sur- 
veys of Orion-KL report percentages of unidentified lines 
between 8 and 15% (Schilke et al. 1997, 2001; Comito 
et al. 2005). Orion-KL is known to have a number of dif- 
ferent velocity components approximately between 2 and 
9kms _1 (e.g., Genzel & Stutzki 1989), and we cannot de- 
termine at which velocity an unidentified line is before- 
hand. Hence, at 690 GHz the frequency range of potential 
molecular lines between 2 and 9kms~ 1 is 16 MHz, com- 
pared to 8 MHz at 345 GHz {Av = Av x v/c). Further- 
more, going to higher frequencies, astronomical spectral 
line studies are rare and identifications that rely only on 
laboratory work are difficult. Hence, identifying spectral 
lines gets more complicated at higher frequencies. 

Figure 8 presents some of the integrated line emission 
maps of the ULs. Since these lines are on average of lower 
intensity, the signal-to-noise ratio in the maps is worse 
than for the identified lines in Fig. 6. The rest-frequencies 
listed for the ULs in Table 1 (and shown in Fig. 8) corre- 
spond to the lines being set to a vi sr of 5kms _1 , implying 
the potential range of frequencies for each line as discussed 
above. In principle, it should sometimes be possible to 
associate unidentified lines with molecular families (e.g., 
oxygen-, nitrogen-, or sulphur-bearing species) based on 
the spatial distribution of the gas (see, e.g., the previous 
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865 /im Orion line studies by Beuther ct al. 2005). How- 
ever, since the quality of the images is rather poor, this is 
a difficult task for this dataset. Nevertheless, we tried to 
associate some unidentified lines with potential molecular 
transitions. Table 1 lists a few tentative candidate lines 
for the ULs, which, however, are only suggestions and not 
real identifications. 

4. DISCUSSION AND CONCLUSIONS 

The presented line and continuum data of Orion-KL in 
the 440yum/690GHz band show the power of high-spatial- 
resolution studies in the submm wavelength bands. The 
measured continuum flux from source / allows us to dif- 
ferentiate between various previously proposed physical 
models: Source / appears to be a rather "normal" pro- 
tostellar object with a SED fitted by a two-component 
model of proton-electron free-free emission below 100 GHz 
plus a dust component contributing to the flux in the 
submm bands. Furthermore, the source SMAl becomes 
more prominent at higher frequencies and is clearly dis- 
tinguishable from the hot core emission. Since SMAl is 
detected neither at cm nor at infrared wavelengths, it 
may be one of the youngest sources of the evolving clus- 
ter. The only weak detection of the hot core at 440 [im is 
puzzling. Although it might be a real feature of the new 
high-spatial- resolution 440 /im continuum data, we cannot 



entirely exclude that it is caused by the poor uv-sampling 
and resulting imaging problems. Keeping the uncertain- 
ties in mind, we find a lower spectral index toward the hot 
core compared with source / and SMAl. This is consis- 
tent with lower spectral indices found toward Orion-KL in 
lower-spatial-resolution single-dish observations that are 
dominated by the hot core emission. 

The spectral line maps trace a warm gas component at 
the center of the region, mainly confined to the sources 
/, SMAl and the hot core peak position. Tempera- 
ture estimates based on the CH 3 CN(37^: — 3&k) K-ladder 
(K = 0...3) lines indicate a warm gas component in the 
regime of 600 ± 200 K. The number of unidentified lines in 
the given setup is large. Potential reasons are discussed, 
of which the two main ones are likely the large potential 
spread in velocity, and thus frequency, of the ULs, and the 
less explored submm wavelength band. 

We would like to thank Peter Schilkc for providing the 
XCLASS software to model the CH 3 CN spectra. Thanks 
very much also to Darek Lis for discussing spectral indices 
in Orion-KL. We also appreciate the helpful comments 
from the referee. H.B. acknowledges financial support by 
the Emmy-Noether-Program of the Deutsche Forschungs- 
gemeinschaft (DFG, grant BE2578/1). 



REFERENCES 



Allen, D. A. & Burton, M. G. 1993, Nature, 363, 54 

Beuther, H., Zhang, Q., Greenhill, L. J., et al. 2004, ApJ, 616, L31 

— . 2005, ApJ in press 

Blake, G. A., Mundy, L. G., Carlstrom, J. E., ct al. 1996, ApJ, 472, 
L49 

Chernin, L. M. & Wright, M. C. H. 1996, ApJ, 467, 676 
Chrysostomou, A., Burton, M. G., Axon, D. J., et al. 1997, MNRAS, 
289, 605 

Comito, C, Schilke, P., Phillips, T. G., et al. 2005, ApJS, 156, 127 
de Vicente, P., Martin-Pintado, J., Neri, R., & Rodriguez-Franco, A. 

2002, ApJ, 574, L163 
Dougados, C, Lena, P., Ridgway, S. T., Christou, J. C, & Probst, 

R. G. 1993, ApJ, 406, 112 
Friedel, D. N., Snyder, L. E., Turner, B. E., & Remijan, A. 2004, 

ApJ, 600, 234 
Genzel, R. & Stutzki, J. 1989, ARA&A, 27, 41 

Gezari, D. Y., Backman, D. E., & Werner, M. W. 1998, ApJ, 509, 
283 

Greenhill, L. J., Gezari, D. Y., Danchi, W. C, et al. 2004, ApJ, 605, 
L57 

Gurwell, M. A. 2004, ApJ, 616, L7 

Gurwell, M. A. & Muhleman, D. O. 1995, Icarus, 117, 375 
— . 2000, Icarus, 145, 653 

Harris, A. I., Avery, L. W., Schuster, K.-F., Tacconi, L. J., & Genzel, 

R. 1995, ApJ, 446, L85 
Ho, P. T. P., Moran, J. M., & Lo, K. Y. 2004, ApJ, 616, LI 
Johnstone, D. & Bally, J. 1999, ApJ, 510, L49 
Lis, D. C, Serabyn, E., Keene, J., et al. 1998, ApJ, 509, 299 
Liu, S., Girart, J. M., Remijan, A., & Snyder, L. E. 2002, ApJ, 576, 

255 



Miiller, H. S. P., Thorwirth, S., Roth, D. A., & Winnewisser, G. 

2001, A&A, 370, L49 
Masson, C. R. 1994, in Astronomical Society of the Pacific 

Conference Series, 87 
Menten, K. M. & Reid, M. J. 1995, ApJ, 445, L157 
Plambeck, R. L., Wright, M. C. H., Mundy, L. G., & Looney, L. W. 

1995, ApJ, 455, L189 
Poyntcr, R. L. & Pickett, H. M. 1985, Appl. Opt., 24, 2235 
Sault, R. J., Teuben, P. J., & Wright, M. C. H. 1995, in ASP Conf. 

Ser. 77: Astronomical Data Analysis Software and Systems IV, 433 
Schilke, P., Benford, D. J., Hunter, T. R., Lis, D. C, & Phillips, 

T. G. 2001, ApJS, 132, 281 
Schilke, P., Groesbeck, T. D., Blake, G. A., & Phillips, T. G. 1997, 

ApJS, 108, 301 

Schultz, A. S. B., Colgan, S. W. J., Erickson, E. F., ct al. 1999, ApJ, 
511, 282 

Scoville, N. Z., Carlstrom, J. E., Chandler, C. J., ct al. 1993, PASP, 
105, 1482 

Shuping, R. Y., Morris, M., & Bally, J. 2004, AJ, 128, 363 
Stolovy, S. R., Burton, M. G., Erickson, E. F., et al. 1998, ApJ, 492, 
L151 

Sutton, E. C, Blake, G. A., Genzel, R., Masson, C. R., & Phillips, 

T. G. 1986, ApJ, 311, 921 
Wilncr, D. J., Wright, M. C. H., & Plambeck, R. L. 1994, ApJ, 422, 

642 

Wilson, T. L., Gaume, R. A., Gensheimcr, P., & Johnston, K. J. 

2000, ApJ, 538, 665 
Wright, M. C. H., Plambeck, R. L., & Wilner, D. J. 1996, ApJ, 469, 

216 



SMA 440^m/690GHz observations of Orion-KL 



7 



[Jy] 



50 






6.8 10 5 


6.805 10 5 .6.81, 10 5 
Rest Frequency (MHz) 


6.815 10 5 




[Jy] 



100 - 



6.9 10 



.905 10 ,6.91 10" 

Rest Frequency (MHz) 



6.915 10" 



Fig. 1. — Vector-averaged lower and upper sideband spectra in the uv-domain on a baseline of ~ 25 m. The data are averaged over an 
hour angle range from —2.6 to +4.5 hours with a total on-source integration time of 176 minutes. UL marks unidentified lines and "?" marks 
tentatively identified lines. The lines at the bottom of the USB spectrum mark the apparently line-free part of the spectrum, which we used 
to produce the 440 /im continuum image (Figs. 3 & 4). 
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Fig. 2. — The top panel shows the uv-coverage of the observations. The bottom panel presents the resulting dirty beam, contour levels 
cover the range ±1, in steps of ±0.125. 



Fig. 3. — Submillimeter continuum images of Orion-KL. Positive emission is always shown in gray scale with contours, negative features - 
due to missing short spacings and the calibration/imaging difficulties discussed in §2 — are presented in dashed contours. The left panel shows 
the previously published 348 GHz continuum data (Beuther ct al. 2004). The contours start at the 2<r level of 70mjybeam and continue 
in 2a steps. The main sources in the region are labeled by name. The middle panel shows an image of the same region, with a limited range 
of uv-data between 30 and 160 kA. The contouring is the same as in the left panel, and the stars mark the four sources identified in the left 
panel. The right panel shows the new 690 GHz continuum data covering a uv-range between 30 and 160 kA again (blanking the only bin 
below 30 kA at 28 kA). The contouring is done in 10 to 90% steps from the peak emission of 10.3 Jy, and the stars again mark the sources 
identified in the left panel. The synthesized beams are presented at the bottom-left of the panels. The different synthesized beams in the 
middle and right panel result from a different uv-sampling, although we applied the same uv-limits. 



Fig. 4. — Original image and model-fits to the 690 GHz continuum data. 1.) Original image, 2.) 2 point source model, 3.) 2 point source 
+ Gaussian model, and 4.) 3 point source + Gaussian model. The contouring is done in ±10 to ±90% steps from the peak emission of the 
three images. The peak values are 10.3, 7.7, 10.2 and 10.6 Jy for panels 1 to 4, respectively. The stars again mark the sources identified in 
the left panel of Figure 3. The synthesized beams are presented at the bottom-left of each panel. 
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Fig. 5. — The SED of source /. The measured fluxes are labeled as filled circles with error bars, and the various lines show different attempts 
to fit the data as labeled within the plot. The 15, 22 and 86 GHz fluxes are taken from Plambeck et al. (1995) and references therein, the 8 
and 43 GHz measurements are more recent values (consistent with the error bars of the previous measurements by Menten & Reid (1995), 
M.J. Reid et al. in prep.), the 345 GHz data point is from Beuther et al. (2004), and the 690 GHz point is from this work. 



Fig. 6. — Integrated line emission maps after continuum subtraction: The top-left panel shows for comparison again the continuum image 
as presented in Fig. 3, and the other 5 panels show representative molecular line maps as labeled within each panel. Positive emission is 
presented in gray scale and full contours, and negative emission - due to missing short spacings - is shown in dashed contours. The contours 
always range from 10 to 90% of the peak emission in each image. The synthesized beams are shown at the bottom-left of each panel. 
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Fig. 7. — CHaCN(37x — 36k) spectra taken toward the CH3CN peak positions associated with source 7, SMA1 and the hot core (see 
Fig. 6). The full lines show the data, and the dotted lines LTE model spectra (adopted temperatures are given in each panel). The data are 
not very sensitive to the temperatures and can be modeled with T varying approximately ±200 K. 



Fig. 8. — Integrated line emission maps of some unidentified lines after continuum subtraction: positive emission is presented in gray scale 
and full contours, and negative emission — due to missing short spacings — is shown in dashed contours. The contours always range from 10 
to 90% of the peak emission in each image. The rcst-frcqucncy for each line set to a v\ SI = 5kms _1 is shown at the top-left of each map. 
The symbols are the same as in the previous Figures, and the synthesized beams are shown at the bottom-left of each panel. 
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a Thc frequency uncertainties are below 1 MHz, mostly even below 0.1 MHz. The frequencies for the ULs are 
derived via setting the v\ aI to 5 kms -1 . 

b FWHM of Gaussian fits to the vector-averaged spectra of the short 25 m baseline (see Fig. 1). 

c Tentative molecular transition identifications which are still very uncertain! 

^UL: Unidentified line. 

e This line is at the edge of the bandpass and thus only partly detected. Therefore, the frequency is even 
more uncertain. 

Line blend. 

s Duc to missing flux, a line fit does not make sense. 
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